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Abstract 

We start from a local and linear spacetime relation between tlie electromagnetic excitation 
and the field strength. Then we study the generally covariant Fresnel surfaces for light rays and 
light waves. The metric and the connection of spacetime are left unspecified. Accordingly, our 
framework is ideally suited for a search of possible violations of the Lorentz symmetry in the 
photon sector of the extended standard model. We discuss how the skewon part of the constitutive 
tensor, if suitably parametrized, influences the Fresnel surfaces and disturbs the light cones of 
vacuum electrodynamics. Conditions are specified that yield the reduction of the original quartic 
Fresnel surface to the double light cone structure (birefringence) and to the single light cone. 
Qualitatively, the effects of the real skewon field can be compared to those in absorbing material 
media. In contrast, the imaginary skewon field can be interpreted in terms of non-absorbing media 
with natural optical activity and Faraday effects. The astrophysical data on gamma-ray bursts are 
used for deriving an upper limit for the magnitude of the skewon field. 
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INTRODUCTION 



Classical electrodynamics is fundamentally linked to the concepts of relativity and space- 
time symmetry. In particular, many high precision tests of relativity are based on elec- 
tromagnetic phenomena. At present, there is a growing interest (both theoretically and 
experimentally) in the search for possible violations of the Lorentz symmetry. The latter 
may naturally occur, for example, in the framework of unified theories of physical interac- 
tions at the Planck length scale, and may then manifest themselves at the scales of usual 
high-energy particle physics in the form of corrections to the established quantum field- 
theoretical models. Even though presently the limits of the standard model of elementary 
particle physics (SM) become visible, see Sec. 10.6 of the new Particle Data Report it still 
remains the basis of our understanding of the strong and the electroweak interactions. At 
the foundations of the SM lays special relativity with the Poincare symmetry group and, as- 
sociated with it, the light cone structure. Although this rigid light cone structure of special 
relativity is assumed as framework for the SM, we know from general relativity theory (GR) 
and from experiments that the light cone structure is in fact "flexible". In other words, in 
GR the metric becomes a space and time dependent field. This alone demonstrates that the 
SM cannot be strictly correct. 

Numerous attempts have been made to take care of this state of affairs and to build up a 
more general framework for the SM and perhaps even to include GR. One of the promising 
approaches is the standard model extension (SME) of Kostelecky and collaborators, see 
I2I; Q]; in which by means of a set of parameters certain violations of the SM are put in 
a quantitative form amenable to experimental tests. The photon sector of the SME has 
been particularly closely studied in [5], where the propagation of electromagnetic waves in 
vacuum was investigated with Lorentz violating terms added in the effective constitutive 
tensor. 

In this paper, we continue the study of electrodynamical phenomena in the framework of 
the so-called premetric approach, i.e., without assuming a special form of the metric structure 
on a spacetime manifold. Such an approach provides a general technique for an appropriate 
physical discussion of the experimental tests of possible violations of the Lorentz symmetry 
and of parity, since at the outset no assumptions were allowed about geometric structures, 
such as the metric or the connection of spacetime. On the contrary, the analysis of the 
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dynamics of the electromagnetic field in this approach yields, via the light propagation, the 
information about the underlying metrical structure of spacetime. Accordingly, this scheme 
can be viewed as a general test theory of the searches for possible deviations from the Lorentz 



symmetry in the photon sector of SM 



approach of Kostelecky and Mewes [s] can be naturally embedded into the framework of 
general premetric electrodynamics. 

Electromagnetic wave propagation is a very important physical phenomenon in classical 
field theory. In general, the geometrical structure of spacetime as well as the intrinsic 
properties and the motion of material media can affect the light propagation. Hence, a 
theoretical analysis of the latter ultimately results in establishing the properties of the 
(genuine or effective "optical" ) metric structure on the manifold. 



It is worthwhile to note that the multi-parameter 



In the generally covariant premetric approach to electrodynamics ll|, |l2|, [is], Q], the ax- 



ioms of electric charge and of magnetic flux conservation manifest themselves in the Maxwell 
equations for the excitation H = {V, Ti.) and the field strength F = {E, B): 

dH = J, dF = 0. (1) 

These equations should be supplemented by a constitutive law H = H{F). The latter rela- 
tion contains the crucial information about the underlying physical continuum (i.e., about 
spacetime and/or about the material medium). Mathematically, this constitutive law arises 
either from a suitable phenomenological theory of a medium or from the electromagnetic 
field Lagrangian. It can be a nonlinear or even nonlocal relation between the electromagnetic 
excitation and the field strength. 

Earlier, we have investigated the propagation of waves in the most general local and 
linear theory. One can conveniently split the constitutive relation into the three irreducible 
pieces which are known as the principal part, the skewon part and the axion part. It is 
important to stress that such a decomposition is intrinsically metricfree, and thus it suits 
nicely for the purpose of the study of the violations of the Lorentz symmetry in the photon 
sector. More specifically, we investigated the influence of the axion and the skewon fields 



on the wave propagation by applying a linear response formalism, see ^|. Subsequently, 
Itin [7 1 studied birefringence caused by the axion field in the Caroll-Field-Jackiw model 



of electrodynamics 



More recently, we investigated the 



waves and its consequences for the metric of spacetime 



airefringence of electromagnetic 
qI. The possible limits of the 



Lorentz symmetry can be seen in a particularly transparent way, both qualitatively and 
quantitatively, including also the case of a possible violation of the charge conservation law 



io|. 



The axion field and its interference with electrodynamics is a well-studied subject, and 
the same applies to the influence of the general principal part of the constitutive tensor 
which was thoroughly analyzed in Sj. However, this is not so in the case of the skewon 
field, a new structure which allows us to exhaust all 36 parameters of a linear spacetime 
relation (exact definitions are given below). The only relevant earlier results are those 



of Nieves and Pal 
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161 ] which were confined to the case of the spatially isotropic and 
constant scalar skewon field. In contrast, we present here some results which hold true for all 
electro dynamical models with an arbitrary local and linear spacetime relation that contains 
a nontrivial skewon part. We use the expression "spacetime relation", since it applies to 
spacetime ( "the vacuum" ) itself, and in order to distinguish it from the constitutive law of a 
material medium. From a mathematical point of view, the general local and linear spacetime 
relation is similar to the constitutive relation for matter with sufficiently complicated electric 
and magnetic properties. Quite naturally, the study and the physical interpretation of the 
wave propagation in a manifold with a general spacetime relation appears to be very close 
to the analysis of the corresponding wave effects in crystal optics. The latter subject was 



extensively studied in the literature, see 
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20|], for example. In particular, our 



analysis reveals the analogy of the effects of the real skewon field to the propagation of waves 
in absorbing media with circular dichroism, whereas the imaginary skewon field appears to 
be analogous to the optical activity tensor. 

If local coordinates are given, with i,j, ... = 0, 1, 2, 3, we can decompose the excitation 
and field strength 2-forms into their components according to 

H = ^ Hij dx' Adx\ F = i Fi^ dx' A dx^ . (2) 

GENERAL LOCAL AND LINEAR CONSTITUTIVE RELATION 

We confine ourselves to the case in which the electromagnetic excitation and the field 
strength are related by the local and linear constitutive law, 

Hij = - Kij''^ Fki . (3) 



The constitutive tensor k has 36 independent components. One can decompose this object 
into its irreducible pieces. Obviously, contraction is the only tool for such a decomposition. 
Following Post |2]|, we can define the contracted tensor of type [{] 



■■= Kii''^ , (4) 
with 16 independent components. The second contraction yields the pseudo-scalar function 

K := Kfc^ = Kki^^ . (5) 

The traceless piece 

^^'■.= K,^-\K6f (6) 

has 15 independent components. These pieces can now be subtracted out from the original 
constitutive tensor. Then, 

^.kl _ (l) . kl . (2) kl . (3) kl /yx 

= + 2 +i/.5f4. (8) 

By construction, •^^^/tjj'^' is the totally traceless part of the constitutive tensor: 

= 0. (9) 

Thus, we split k according to 36 = 20 + 15 + 1, and the [2] tensor (^^Kj^'^' is subject to the 
16 constraints ([9]) and carries 20 = 36 — 16 components. 



One may call ^^^Kij'''' the principal, or the metric-dilaton part o: 
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the constitutive law. 
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2J|. We further 



Without such a term, electromagnetic waves are ruled out, see 
identify the two other irreducible parts with a skewon and an axion field, respectively. 
Conventionally, the skewon and the axion fields are introduced by 

^i^ = -\t^^, c^ = ^'^- (10) 

The standard Maxwell-Lorentz electrodynamics arises when both skewon and axion van- 
ish, whereas 

= Aor/./. (11) 



Here Aq = ^Jso/l^o is the vacuum impedance. A spacetime metric gij is assumed on the 
manifold, and with g := detgij one defines rjijki := yj—ghiki and r/jj^' = rjij^n 9^'' Q""^ ■ It has 



been shown recently [9j that taking the hnear spacetime relation for granted, one ends up 



25|. 



at a Riemannian lightcone provided one forbids birefringence in vacuum, see also 

Along with the original K-tensor, it is convenient to introduce an alternative representa- 
tion of the constitutive tensor: 

x'''''-=\e''"'''ti^n'''. (12) 
Substituting ^ into ffT^ . we find the corresponding decomposition 

with the principal, skewon, and axion pieces defined by 

(2) ijkl _ \ ijmn (2) kl _ _ ijm[k i l\ /i c\ 

A ~ 2 "'mn — fc ; K'-^ I 

{3)^ijkl ^ ^^ijmn = }_^ijkl^_ ^^g^ 



Using the S-identity and the K-identity derived in 23|, we can verify that (^^x is skew- 
symmetric under the exchange of the first and the second index pair, whereas ^^^x is sym- 
metric: 

{2)^ijkl ^ _(2)^klij ^ (l)^ijkl ^{l)^klij_ ^-^j-^ 



Space-time decomposed constitutive relation 



Making a (1 + 3)-decomposition [13| of covariant electrodynamics, we can write H and 
F as column 6-vectors with the components built from the magnetic and electric excitation 
3-vectors Ha,'D'^ and electric and magnetic field strengths Ea,B°', respectively. Then the 
linear spacetime relation ([3]) reads: 





Here the constitutive tensor is conveniently represented by the 6-matrix 



K 




X 



IK 




(18) 



(19) 



The constitutive 3x3 matrices A, B, C, V are constructed from the components of the original 
constitutive tensor as 

^ta ^OaOt ^ 1 ^^^^ ^^^^ ^cdef ^ (2O) 

C\ := \ hcd X'""'"' , -Da' ■■= I eacd X''"" ■ (21) 
If we resolve with respect to x, we find the inverse formulas 

^OaOb Aba ^nbcd ^abe ^cdf ^00^ 

X — X — e e ' Dfe , [ZZ) 

^Oabc ^ ^bcd ^ ^abOc ^ ^abd qc^ _ ^23) 

The contributions of the principal, the skewon, and the axion parts to the above consti- 
tutive 3-matrices can be written explicitly as 

jab ^ _^ab _^abc^Q^ (24) 
Bab = l^ab + ^abc , (25) 

C\ = i\ - - S^b ^c") + a 51 (26) 

Va' = 7'a +(»^a'-5^^/)+«e (27) 

The set of the symmetric matrices e"'' = e'^"' and fi~,^ = fi^^ together with the traceless matrix 
(i.e., 7^c = 0) comprise the principal part ^'^^■y^^i^^ of the constitutive tensor. Usually, 
e"'^ is called permittivity tensor and reciprocal permeability tensor ("impermeability" 
tensor), since they describe the polarization and the magnetization of a medium, respectively. 
The magnetoelectric cross-term is related to the Fresnel-Fizeau effects. The skewon 
contributions in fl2^ and fl25l) are responsible for the electric and magnetic Faraday effects, 
respectively, whereas skewon terms in 0261) and (1271) describe optical activity. 

GENERAL FRESNEL EQUATIONS: WAVE AND RAY SURFACES 



22,M,M. In the Hadamard 



Here we briefly summarize the results of previous work 
approach, one studies the propagation of a discontinuity in the first derivative of the elec- 
tromagnetic field. The basic notions are then the fields of the wave covector and the ray 
vector that encode the information about the propagation of a wave in a spacetime with a 
general constitutive relation. 
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Wave surface 



The crucial observation about the surface of discontinuity 5* (defined locally by a function 
$ such that $ = const on 5*) is that across 5* the geometric Hadamard conditions are satisfied 
for the components of the electromagnetic field and their derivatives: [Fij] = 0, [diFj^] = 
qi fjk, [Hij] = 0, [diHjk] = qihjk- Here qi := di^^ is the wave covector. Then using the 
Maxwell equations ([1]) and the constitutive law ([3]), we find a system of algebraic equations 
for the jump functions: 

x'^''g,/fc/ = 0, e^^-^'g,/,, = 0. (28) 

Solving the last equation in fl28l) by means of fij = qiaj — qjQi, we finally reduce fl28l) i to 
^ijki q.q^d^ = 0. This algebraic system has a nontrivial solution for only if the determinant 
of the matrix on the left hand side vanishes. The latter gives rise to our generalized covariant 
Fresnel equation 

G''''\x)mjqkqi = 0, (29) 
with the fourth order Tamm-Rubilar (TR) tensor density of weight +1 defined by 

G'^'^ix) ■■= ^ e^np, erstu X^'^'^' X^''"^" X'^'^*" • (30) 

It is totally symmetric, ^*-''^'(x) = Q'^''^^^^ [x) ■, and thus has 35 independent components. 
Different irreducible parts of the constitutive tensor ([p contribute differently to the 



general Fresnel equation. A straightforward analysis [23|, |2J| shows that the axion piece 
drops out completely from the TR-tensor, whereas the two remaining irreducible parts of 
the constitutive tensor contributes to (|5U]) as follows: 

Q^3^\^) = g^^f'^Wx) + (^^x™^''"'-'' ^J: ^r!^ ■ (31) 



Ray surface 



A ray can be defined by a vector field s = s^di that is dual to the wave covector q = qi dx 



27 



28|: 



and to the wave front in the following sense 

sj/i = 0, sJ/ = 0, sjg = 0. 



(32) 



Here the 2-forms h = ^ hijdx^ A dx^ and / = ^ fijdx^ A dx^ represent the jumps of the first 
derivatives of the electromagnetic field. Equations (|32|) are metric-free. The system of the 
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first two equations can be analysed in complete analogy to the system (128|) . As a result, one 
arrives at the "dual" Fresnel equation imposed on the components of the ray vector: 

Gijki s'sh's' = 0. (33) 

Like the TR-tensor density fl3UI) . the totally symmetric tensor density of weight —1 is con- 
structed in terms of the components of the constitutive tensor: 

A . ^ mnpq rstu O C C 

yijkiyX) ^ Xmnr{iXj\ps\kXl)qtu 

4 4! '^mnr(i ^j\pvw A ^ghs\k ^l)qtu A. ■ V / 

As before, a complete symmetrization has to be performed over the four indices i, j, I, with 
the vertical lines separating out those indices that are excluded from the symmetrization. 
Here we introduced a "double dual" of the constitutive tensor density as 

Xijkl T ^ijmn ^klpq X 

(35) 

The contraction of this new object with the original constitutive tensor yields 

1 ..1 

2 Xklmn X ~ 2 "'mn • J 

Using (fT3i) . we find the corresponding decomposition 

Xijkl = ^^^Xijkl — '^hj'm[k ^l]^ + hjkl (37) 

Substituting this into fl34l) . we arrive at the skewon contribution to the ray surface: 

G^MX) = G^Jkl{^'h) + ^'^Xmmj h"^ h)" ■ (38) 

SKEWON FIELD: DIFFERENT PARAMETRIZATIONS AND WAVE PROPAGA- 
TION EFFECTS IN VACUUM 

The main aim of this paper is to investigate the possible infiuence of the skewon field 
on the electromagnetic wave propagation. For simplicity, we assume that the principal part 
of the constitutive tensor is determined by the Maxwell-Lorentz law ( JTTl) . that is. 



= Ao (f'g'' - g^-'g^') , ^^^Xn.^n, = - ^ ^9.^ng,i ' g^.g^) . (39) 

However, the skewon structure will be kept as general as possible. As a preliminary step, 
we note that direct computations yield 
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'Double vector' parametrization 



The skewon, as a traceless tensor field of type [\], has 15 independent components. Given 
the spacetimes metric gij, we can decompose the skewon field into an antisymmetric and 
a symmetric part. The latter can be conveniently constructed from a pair of arbitrary 
(co)vector fields f * and Wi , so that we find eventually 



'i^ = + WiV^ - - 51 (wv). 



(41) 



Here aij := ai^Qkj = —aji is an arbitrary antisymmetric tensor and (wv) = WiV^. Such a 
parametrization produces almost a general skewon field, since the number of independent 
components is here 14 = 6(ajj) + 4(f*) + 4(wj). 

Substituting ( |4T|) into ( !3T1) and ( !38l) . and taking into account ( !39l) . we find 



(1) 



m CP n 



Xm{i\n\j pk pi) 



, (42) 
• (43) 



Here 6* := a/, 6j := Oj-' f j and w"^ = WiW\ v"^ = ViV\ 

Adding the principal terms fHOj) to the contributions of the skewon fH2|) and fH3|) . we 
see that neither the wave surface nor the ray surface decompose into the product of cones 
in general. In other words, birefringence is absent. The last terms in (1121) and in fH5]) are 
responsible for that. Let us investigate those particular cases of the skewon field in the 
'double vector' representation pTl) which do admit the birefringence phenomenon. 



Purely antisymmetric skewon 

When the skewon field is represented by its antisymmetric part only, i.e.,^i^ = (which 
is obtained by putting both = Wi = 0), then the analysis of (131!) shows that the light 
propagation is birefringent with one light cone defined by the spacetime metric g^^ and the 
second "optical" metric 

4~>=9''-^a^^aJ. (44) 

This result actually holds true for any other parametrizations of the skewon field which will 
be considered later. At the same time, inspection of fl38l) yields the decomposition of the 
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ray surface into the light cone of the spacetime metric gij and a second cone defined by 

(2) 1 

9ij 9ij Tf '^mi'^ j- {^^) 

Ahhough fj44|) and fj45l) look very similar, we notice that they are not mutual inverse, and 
hence, they do not determine a unique geometric structure on spacetime. 



Purely symmetric 'vector' skewon 



Assume now that the skew-symmetric skewon part is absent, which is evidently achieved 
if aij = V[iWj]. Then 6* = {l/2)[w\vw) — v^w"^], hi = {l/2)[vi{vw) — Wiv'^], and we find the 
TR-tensor in the birefringent form 



— (ij kl) 

-g g+g-. 



Here the two light cones are defined by the two optical metrics 



{vw) ± J{vwY - 4A§ 



(46) 



(47) 



Here {yw) = Viw\ Obviously, in order to have the wave propagation along separate light 
cones, the skewon field must satisfy the condition 



{vwf > AXl 



(48) 



When this condition is violated, the Fresnel wave surface is of 4th order and birefringence 
is absent. For the equality sign in fHHj) . the birefringence disappears and the two light cones 
coincide. However, one can straightforwardly prove that for [vwY = 4Aq the optical metric 
has an Euclidean signature. Hence there is no wave propagation in this case. The type of 
the signature correlates with the sign of the determinant of the metric: the latter is negative 
for the Lorentzian case and positive for the Euclidean one. 
A direct computation yields 

det [g'^ - kv^'w^^) = (det g'^) \l-k {vw) + ^ [{vwf - v^w^] \ . (49) 



Comparing this with fj47|) . we identify the coefficient as k± = (vw) ± (vw)'^ — 4Aq /2Aq 
and finally obtain 



(det g'i) = (det g'^) ^ [{vw)^ - vW - 4X1 



(50) 
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Since the spacetime metric is Lorentzian, its determinant is negative, {det g^^) < 0. As a 
result, we conclude that the optical metrics will be also Lorentzian only if 



[vw 



2 2 ^ ^ \ 2 
V W > 4Ao. 



(51) 



Note that the two conditions (!5T]) and fHHj) look pretty similar. In fact, we can verify that 
only one of them is an essential inequality. Indeed, when the two vectors are either both 
timelike {v"^ > 0,w^ > 0) or both spacelike {v"^ < 0,w^ < 0), the new condition (1511) makes 
the constraint (148|) redundant, whereas when at least one of the vectors is null, (ISTIl and 
fHHj) coincide. Finally, when one of the vectors is timelike and another spacelike, the only 
essential condition is (jUj) and (ISTl) is redundant. 

With these conditions, we can classify the effects of the skewon field on wave propagation. 
The spatially isotropic skewon field of Nieves & Pal 15|, llGj, for example, belongs to the 
subcase when the two vectors w and v are proportional to each other. Indeed, let us consider 
the (fiat) Minkowski spacetime with the Lorentzian metric Qij = diag(c^, —1, —1, —1), and 
let us put = — f ' 



6Qy2S/c with some scalar function S of dimension [Aq]- Then, from 



pTj) . we recover the isotropic skewon field 



S_ 
2 



V 



-3 
10 
10 
1 



(52) 



As a result, the spacetime relation becomes 

V = 60 5"^ El, + (-5 + a) 5^ 



(53) 
(54) 



Accordingly, in the special case when skewon and axion become constant fields, one can 
speak of the 4 electromagnetic constants for a vacuum spacetime with spatial isotropy: The 
electric constant Eq, the magnetic constant /xq, the spatially isotropic part 5* of the skewon 
, and the axion a. 

In this special case we evidently find {vw)'^—v'^w'^ = 0. Hence fl5U]) yields that both optical 
metrics have a wrong Euclidean signature. This is in complete agreement with our earlier 
observations (23! in which the absence of wave propagation was reported on the background 
of the isotropic skewon field of Nieves & Pal. 
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However, now we obviously obtain a far more general result: whereas Nieves & Pal 
have considered only a very special configuration with the spacelike vectors v and w being 
proportional to each other, the wave propagation is actually also absent for any ~ f * and 
any Riemannian spacetime metric Qij. Then again {vwY — v'^w'^ = and hence (detg^^) = 
-{detg'^)Xlkl > 0. 



General parametrization 

It is impossible to construct a general parametrization of the skewon field by using covari- 
ant 4-dimensional objects. Instead, we have to switch to the space and time decomposition 
techniques and to allow for the use of 3-dimensional (3D) tensor objects. Then we can write 
the required parametrization as follows: 



(55) 



V 



Here and rib are 3D (co)vector fields, whereas Sb"' is a 3D tensor of type [\]. From now on, 
the indices from the beginning of the Latin alphabet denote the spatial 3D components, i.e., 
a, 6, c, . . . = 1, 2, 3. The count of the independent components, namely 3 + 3 + 9 = 15, shows 
that we indeed have the most general representation of the skewon field. It is straightforward 
to establish the dimensions of the different pieces of skewon field: We find that [sb"'] = [Aq], 
[ub] = [eo] and [m"] = l/[/io]. 

With this general parametrization, the linear spacetime relation (llSp can be written 
symbolically as follows: 

f n-/ \ /\ / _ \ 

(56) 





Here the linear operator on the right-hand side explicitly reads: 



b„ c 



^abc 



-eabcfn" 

Sb" - 5fe 



/ ,b 



a 



K 

6^ 



(57) 



Now we again specialize to the Maxwell-Lorentz principal part fl5Ul) . Moreover, since we 
are specifically interested in the skewon effects, we will assume for concreteness that the 
spacetime metric is the flat Minkowskian one with Lorentzian signature, Qij = (c^, — 6ab) = 
diag(c^, — 1, — 1, — 1). The Euclidean 3-metric tensor 6ab should not be confused with the 



13 



Kronecker delta: The components of the latter are always either or 1, whereas the former 
has these values in Cartesian coordinates only. Then one can further decompose the 3D 
skewon tensor into its symmetric and skewsymmetric parts: 

st^ = Uh'' + e^Zc. (58) 

Here Uab = Uba, and the 3D indices are raised and lowered by means of the Euclidean 
3D metric, i.e., Ub"" =: UbcS"''^ and e^'^'^ := ebde^'^'^S'^^- The covector field Zc describes the 
antisymmetric part of the 3D tensor skewon. Obviously, both irreducible parts of the tensor 
skewon have the same dimension [ub"'] = [zc] = [Xq]. 



Space-time decomposition and wave and ray surfaces 



The analysis of the wave propagation in the presence of the general skewon field fl55|) 
is not possible in the above covariant 4-dimensional framework. Since the space and time 
components are explicitly separated in the general parametrization ([S5D, we need a corre- 
sponding approach to the wave and ray surfaces in which a space and time decomposition 
has been performed. 

Namely, let us denote the independent components of the TR-tensor fl30l) as follows: 



M 



M' 



ab 



M' 



abc 



M' 



abed 



--0000 



det^. 



4 = -ebcd {a''" A"" C^, + A"^ A'' V, 



d\ 
e j 5 



OOab ^ _ j^iab) Y^^d^Y ^ ^^ey _ ^Q'^ ^ + P/)(C'^, + V, 

+{C'^c + V,^){A'^''C^^ d + Vd^^A^^') - C^dA'^X^K 



-V^^'^A^^^Vd'^ - A'^'C^WJ'^ + (A^'^^^A'^' - ^^("^^)") Bde 

^gOabe ^ ^de{e\ \j3^^^j[c^b) _ 1)^"^^)/ ) 

+BfdiA'''^ C{ - A^^X'l) + 1)/) Vj + P/ C'l 

r 1 

gabcd _ ^ef{c^\gh\d 



- 4"^) H — C°- T)^ 



Then, in decomposed form, the generalized Fresnel equation (1291) reads 



(59) 
(60) 



(61) 

(62) 
(63) 



q^M + qtqa + q'^q^qb M^' + q^q^qbqe M"''^ + q^qbqeqd M''""' = . 



abe 



rabed 



(64) 



Similarly to (!59|l - (l63|) . we can decompose the 4-dimensional tensor density (IMl) into a set 
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FIG. 1: Fresnel surface for a real skewon of the electric Faraday type. It has the form of a 
toroid (depicted with a cut). Here Ua = nd^ with = 10; we use the dimensionless variables 
X := cqi/qo, y := cq2/qo, z := cqs/qo. 



of 3-dimensional objects: 



M 



ab 



^0000 = det i3 , 

4 GoOOa = —e'"^'^ {Bba Bee T^d + ^ab ^ec d) i 

6 Gmab = \ Biab) [iC^y + {Ve'f " {C^ d + ^^Z) {C\ + V^") 
HCd + Vf){Be(aVb)'' + C\aBb)c) ' V/B^iaVb)' 
—Ci^aBbyC^d — BdcT^{aC'^b) + {B(ab)Bdc — Bd(aBb)c'j -4' 



dc 



Mabc '■= ^Goabc = ^de{c •A.'^^ {B ab) f — CaB, 



+Af\Bab) Vf'- Blf\aVbf) + Vj C\) + C^a Vbf V 



f 



-Bab)A^'-Va'C\-^ 



(65) 
(66) 



(67) 

(68) 
(69) 



Mabcd '■— Gabcd — ^ef{c^\gh\d-A!^^ 

Then the Fresnel type equation for the quartic ray surface (l33ll in decomposed form reads 



M + S'' Ma + S'^S' Uab + S'^S'^S'S' Mabc + s' Mabcd = . (70) 

Now we are in a position to study the effects induced by the different irreducible parts of 
the skewon field. 



Skewonic electric Faraday effect 

When only the part of the skewon is nontrivial, whereas = and Sf,° = 0, the 
electric constitutive relation reads = EqE"- + e'^^'^nbE^. We will call this the skewonic 
electric Faraday effect, since it describes the electric Faraday effect in optically active media 
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FIG. 2: Fresnel surface for a skewon of the electric Faraday type. It has the form of a spheroid 
for the large purely imaginary skewon n. Here = n5^ with nP = —10; we use the dimensionless 
variables x := cqi/qo, y := 0^2/90, z := cqs/qo- 



for purely imaginary 12|, |20|, ISjJ. Then, from ( I59l) -( l63l) . we find the coefficients of the 
Fresnel equation as M° = 0, M"*'^ = 0, whereas the nontrivial components read 



M 



- £0 {^l + 



n 



M' 



abed 



ah 



M' 



Here, as usual, = 5°'^ni, and = n°-na. 



Bo 



2\l 6"'' + (5' 



(71) 
(72) 



On the other hand, according to flB51) - fl^^ . the analysis of the ray surface yields Ma = 0, 



abc 



0, and 



M = -, 



M, 



ah 



Oah, 

fJ'O 



M, 



abed 



(73) 



As a result, we obtain the following Fresnel surfaces for the wave covectors and the ray 
vectors, respectively: 



{q,/cY (1 + rf) - {q,/cf [2q^ + rfq' - inqf] + {q^ 
{s^cf - 2{socY + s 



+ insf 



0, 
0. 



(74) 
(75) 



Here we denote q^ = qaq"", s"^ = SaS", (nq) = Uaq"', (ns) = UaS"'. Moreover, we use the 
dimensionless skewon vector Ua := na/sQ. 

With the help of the nontrivial vector n, we can split the wave covector qa = qa + i^'to 
its transversal and longitudinal parts, with q^ = qa — n^'^njf'qi, and glj = n~'^nji'qb- Then 
the Fresnel equation fl7^ is rewritten in the form 



i.%lcy - qi (go/c)ll + n') -qi + -2(go/c)^ + 2qi + 



(76) 
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Accordingly, in the plane transversal to the skewon vector Ua-, we find birefringence with 
the Minkowski metric and the second optical metric g'^^ = diag((l + ?T^)/c^, — 1, — 1, — 1). 
However, the Fresnel surface does not decompose into a product of the two light cones. 
Although, for = 0, the wave covector surface is a sphere, it becomes a toroid already for 
a small positive > 0. The latter looks like a very thin and narrow belt for small values 
of n^. However, it becomes a thick toroid with the inner radius decreasing to zero for large 
values of n"^. A typical Fresnel surface for this case is depicted in Fig. [H Physically, the hole 
of this toroid corresponds to the directions in which waves cannot propagate: In spherical 
coordinates, the wave propagation occurs for ^ ^ ^ — ^^Oj where the spherical angle 6 
is counted from the direction of the skewon vector and the limiting angle is determined 
from sin 9o = 2/{l + VTTW). As we see, the hole disappears when — oo. 

This qualitative characterization of the wave propagation can be completed by a numerical 
estimate for the velocity of light. The phase velocity v of the wave propagation is defined by 
the components of the wave covector according to qa = qoka/v, where ka are the components 
of the 3D unit covector. Usually, the mean value of the velocity is of interest that is obtained 
after averaging over the directions of propagation and the polarizations. For the case under 
consideration, a direct calculation yields 

Thus, since is positive, the velocity of the wave propagation turns out to be less than c. 
Making use of the data on the gamma-ray bursts and following [S], we derive as an upper 
limit for the skewonic electric Faraday effect < 3 x 10~^^eg. 

The existence of a hole in the wave covector surface confirms our earlier conclusion about 
the dissipative nature of the skewon. In particular, when we continue to the purely imaginary 
values of n so that the square becomes negative, < 0, the topology of the Fresnel surface 
changes drastically. The hole disappears and the wave covector surface assumes the form of a 
spheroid around the origin. More exactly, there is a single spheroid elongated in the direction 
of the skewon vector Ua for n"^ > —1, and two (embedded one into another) spheroids when 
— 1 < < 0. The average radius of the inner spheroid in the latter case is smaller than 
1 which, in physical terms, means that the velocity of wave propagation is greater than 
c. Although the negative < corresponds to the well known electric Faraday effect in 
crystals, such a superluminal behavior of waves for — 1 < < was never reported in the 
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FIG. 3: Fresnel surface for a skewon of the electric Faraday type. It has the form of two concentric 
spheroids (depicted here with cuts) for a small purely imaginary skewon n. Here ria = n5^ with 
ri^ = —0.5; we use the dimensionless variables x := cqi/qo, y := cq2/qo, z := cq-i/qo. 



literature, to the best of our knowledge. In Fig. [2] and Fig. [3], we depict typical Fresnel wave 
covector surfaces for ff = —10 and = —0.5, respectively. The wave surface for = — 1 
looks qualitatively like Fig. [2] with the top and the bottom of the spheroid at z = ±1. 



Skewonic magnetic Faraday effect 

In the opposite case, when the m° part of the skewon is nontrivial, whereas = and 
Sb"" = 0, the magnetic constitutive relation reads Ha = Ba + eahci^^B'^- We call this case 
the skewonic magnetic Farada y e ffect, since it corresponds to the magnetic Faraday effect 



when is purely imaginary 

UJ, 120, |32| |. From fl59|) - fl63|) . in full analogy to the previous 
subsection, we find that = 0, M"'"^ = 0, whereas the nontrivial components are 

M = -eI, M^^ = 26oXl 5''\ M^^"'^ = - —S^''^ (Xl 5^'^^ + c'^ rrfm'^A . (78) 

/ l.r\ ^ -' 



/Wo 



Analogously, from we find for the ray surface that Ma = 0, Mabc = 0, and 

1 / 1 



M 



M, 



1 



+ m 



Mabcd — £0 Xq 6(^abScd), 



ab 



/Wo 



2Xl6ab + c-' (5, 



abTJl -irLairLb 



(79) 
(80) 



Here again = Sabm^ and = m"-ma. If, similarly to the previous subsection, we 
introduce the dimensionless skewon vector by m"' := m'^/io, we end up with the following 
Fresnel wave and ray surfaces, respectively: 



(go/c)^-2(go/c)^g^ + g^ + {mqf 



0, 



(81) 
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FIG. 4: Fresnel surface for a skewon of the magnetic Faraday type. It has two branches that are 
both hyperboloids for the large purely imaginary skewon m (depicted with a cut into half). Here 
rff' = miJf with nf' = —10; the dimensionless variables are x := cqi/qo, y := cq2/qo, z := cq^/q^. 



{socy (1 + m^) - (soc)2 [2s^ + m^s" - (ms)^] + {s^f = 0. (82) 

Similarly as before, {fnq) = friaq"' and {ms) = rias"'. Repeating the calculation of the velocity 
of the wave propagation for this case, we find the mean velocity Kv"^ >= (? . As a result, 
the analysis of the gamma-ray bursts data does not impose any limit on the pure skewon 
field of the magnetic Faraday type. 

As one can see immediately, the wave covectors are only trivial in the plane transverse 
to the direction of the skewon vector, i.e., when (mg) = 0. The waves propagate in that 
plane without birefringence along the light cone defined by the standard Minkowski metric, 
since then [(go/c)^ — g^]^ = is fulfilled. In contrast, a purely imaginary skewon vector, 
with < 0, yields more interesting results. In this case, the mentioned circle represents a 
particular configuration that turns out to be the intersection curve of the two branches of 
the Fresnel surface. The typical Fresnel wave covector surfaces for the cases = — 10 and 
= —0.5 are depicted in Fig. |4]and Fig. [5], respectively: Both branches are hyperboloids 
or spheroids, intersecting in the z = plane. The wave surface for rff' = —1 is qualitatively 
different: It is represented by paraboloids as shown in Fig. [61 

Comparison of the results of the two last subsections demonstrates that the electric and 
magnetic Faraday skewon effects are dual to each other in the following sense: When we 
interchange the wave and the ray vectors q ^ s and simultaneously the skewon vectors 
n m, then the wave surface (fMl) has the same form as the ray surface ( l82l) and vice 
versa: The ray surface fITS]) has the same form as the wave surface ( 18 ip . In this sense. 
Figs, min] show the ray (wave) surfaces dual to the respective wave (ray) surfaces depicted 
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FIG. 5: Fresnel surface for a skewon of the magnetic Faraday type. It has two branches that 
are both spheroids for a small purely imaginary skewon m (depicted with a cut into half). Here 
m"' = m^f with = —0.5; the dimensionless variables are x := cqi/qo, y := cq2/qo, z := cq^/qo. 



in Figs. [T] and [21 with the mentioned replacement of n ^ m and with the interchange 
between the dimensionless wave covector variables x = cqi/qo, y = cq2/qo, z = cqs/qo and 
the dimensionless ray vector variables x = s^/s^c, y = s'^/s^c, z = s'^/s^c. 



Skewonic magneto- electric optical activity 



Let us now assume that only the tensor part of the skewon is present in (155|) , whereas = 
and Ua = 0. We call this a skewon of the magneto-electric type, since it corresponds to the 



12 



33 



34|. 



case of the natural optical activity in matter for the purely imaginary Sb'^, see 
Furthermore, we consider the effects of the symmetric and skew-symmetric parts fl58|) . In 
the absence of the skew-symmetric tensor part {z^ = 0), we find M"- = 0, M"'^'^ = and 

/Ag + \IU — Ud ]Uc — Ue U' ' ' 



M = -el 



(83) 
(84) 



As a check of the consistency of our formalism, we can verify that the isotropic case uj^ 
— S5a^ reduces to the birefringent case fHBl) and fHTj) with w'^ = — = 5q V2S/ c. 
As to the ray surface, a direct computation yields Ma = 0, Mabc = and 



C2 6aa'Sw M"'"', Mabcd = - Saa'SwScc'Sdd' M^'''^'"' . (85) 



Accordingly, the form of the Fresnel ray and wave surfaces turns out to be exactly the same. 
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FIG. 6: Fresnel surface for a skewon of the magnetic Faraday type. It has two branches that are 



Here = m^f; we use the dimensionless variables x := cqi/qo, y := cq2/qo, z := cq^/qQ. 



A direct calculation of the mean velocity of the light propagation now yields 



As usual, we again introduce the dimensionless skewon variable uj' = Ua'^/Xo. Then, simi- 
larly to the above cases^we find as upper value for the symmetric tensor skewon field from 
gamma-ray burst data [5] uj'uh'' < 7 x 10"^'^ Ag. 

The Fresnel wave covector surface has now quite a different form as compared to the 
two cases above. For concreteness, let us take a tensor skewon with only one nontrivial 
component, = u (^6l ^2 + 6^ 6fj. Then, for small values of the skewon m < 1, the Fresnel 
wave surface is depicted in Fig. [71 As a comment to this figure, let us recall that for a 
vanishing skewon u = we have a pure vacuum spacetime relation, and the Fresnel wave 
surface is then a sphere. With u ^ 0, this sphere degenerates to a pair of highly deformed 
intersecting toroids that, for extremely small values of u, are nearly covering the sphere. 
But with growing u, the toroids becomes thin and thinner. The intermediate situation is 
actually depicted in Fig. [71 Since we do not expect large skewon fields in general, we here 
limit ourselves to the case of a small u. 

In the complementary situation, when m° = 0, = and the symmetric part is absent 
in f[58|) . uj' = 0, we obtain a particular case described in Sec. : Birefringence with the 
Minkowski light cone and the second optical metric 



both paraboloids for the purely imaginary skewon with m 



1 (depicted with a cut into half). 




(86) 




(87) 
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FIG. 7: Fresnel surface for a skewon of the magneto- electric optical activity type. It has two 
intersecting toroidal branches for the real skewon u = 0.8. We use the dimensionless variables 
X := cqi/qo, y := cq2/qo, z := cgs/go- 

Unlike as with the symmetric tensor skewon, here the optical metric is always Lorentzian 
with the determinant (detglf) = — (? j (1 — z^/Aq)^. This means that the waves propagates 
along both light cones (except for a skewon satisfying = 6°'^ZaZh = Aq when the optical 
metric becomes degenerate). A straightforward computation of the mean speed of the wave 
propagation now yields 



and, accordingly, the upper limit for the skewon field from the gamma-ray data is again 



In the previous section we studied the wave propagation on a vacuum spacetime described 
by the principal part W-^/'^i'^'- of the constitutive tensor ( 139|) . We found that the characteristic 
effect of the skewon field in that case was the emergence of holes in the wave covector surface. 
In physical terms this means the complete damping of the wave propagation in certain 
directions. Here we will briefly study the wave propagation in an anisotropic dielectric 
medium and demonstrate that a similar effect occurs in the presence of a skewon field. 

The constitutive relation for anisotropic matter is most easily formulated in terms of 
the 3 X 3-matrix parametrization (fT9l) - (!2Ti) . In the absence of a skewon field, an arbitrary 




(88) 



z^<7x 10-27 A2 



SKEWON EFFECTS IN MATTER 
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dielectric medium is described by the spacetime relation 



-EqE 



ab 



with 5° 



£i 



£2 

V £3 ; 



Bab = Sa 



The matrices Cb = and Va'' = 0. Then, from (I59l)-(l63!), we find M'^ = 0, M 



Tabc 



(89) 



0, 



M 



abed 



^{ab^cd)^2j^^^ and 



ab 



( 



ei{e2 + e-i) 
£2(^1 + £3) 
£3(^1 + £2)/ 



(90) 



As a result, the Fresnel equation for the wave surface (IMl) can be recast into the simple form 



+ 



£2 ^2 



+ 



(?q^ — Ex ql (?q^ — 62 ql c^q^ — £3 ql 



0; 



(91) 



as in the previous section, = qaq"- 

As to the ray surface, we can immediately verify from (!65!) - (!69!) that Ma = and Mate = 0, 
whereas M = fi^^, Mabcd = £o>^lS^ab£~^cd) £i£2£3, and 



/ 



Mab = -^ 



\ 



\ 



(92) 



£2+^3 

61+63 
61+62 J 

Then, from (170|) . one can straightforwardly derive the dual equation of the Fresnel ray surface 
for the components of the ray vector s 



-.1^2 



-.2\2 



.3\2 



+ 



+ 



0. 



(93) 



eis2-c2(s0)2 £2S2-C2(s0)2 £^s^-c^{sOy 

Typical wave and ray surfaces are depicted in Figs. [8] and [9], respectively. Both surfaces 
consist of two non-intersecting branches touching each other exactly in 4 points. 

When the skewon field is present, the mentioned picture becomes more nontrivial. In 
particular, for the case of the isotropic skewon ( l52ll . the constitutive relation is modified by 
adding to fISIJl) the nonvanishing 3x3 matrices 



Ca 
h 



S6l 



Va' = S6l 



(94) 
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FIG. 8: Fresnel wave covector surface for an anisotropic dielectric medium with £i = 39.7, £2 = 
15.4, £3 = 2.3. There are two branches, the outer part of the surface is cut into half in order to 
show the inner branch; we use the dimensionless variables x := cqi/qo, y := cq2/qo, z := cq^/qo. 



Then, a quick computation shows that all the M-coefficients remain the same except that 
fl90|) is replaced with 



5l(52 +£3 -45 ) 












82(61 +es-4S^) 

63{ei+62~AS^) J 



(95) 



Here we introduced the dimensionless skewon variable S := S/Xq. 

Analogously we can verify that the dual M-coefficients are the same as in the skewonless 
case except that ( l92i) is changed to 



M, 



ab 



X2 



^62 + 63- 45' 



V 

















El + 82- 45^ J 



(96) 



Then the equations of wave and ray surfaces cannot be represented any longer in the 
compact form (!9T]) and (p3|) . The skewon affects the form of the surfaces by creating typical 
holes that correspond to the directions of wave covectors along which no wave propagation 
occurs. More specifically, the original two branches of the wave and ray surfaces, depicted 
in Figs. [HlEl merge and form a single non-simply-connected surface. In the process of such 
a merging, the four points where the original two branches touched each other become the 
wormholes through which one can move from the outer surface into the inner one. Such a 
typical skewon effect is shown in Fig. [TOl 



24 



FIG. 9: Fresnel ray vector surface for an anisotropic dielectric medium with £i = 39.7, £2 = 15.4, 
£3 = 2.3. There are two branches, 1/4 part of the outer surface is cut off in order to show the 
second inner branch; we use the dimensionless variables x := s^/s^c, y := s'^/s^c, z := s'^/s^c. 




FIG. 10: Fresnel wave covector surface for an anisotropic dielectric medium with £1 = 39.7, 
£2 = 15.4, £3 = 2.3 in presence of a skewon. The two original branches (cf. with Fig. [^j 
are now merged into a single surface with the wormholes replacing the original four points 
where the branches touched. The surface is cut into half; we use the dimensionless variables 
X := cqi/qo, y := 032/^0, z ■= cqs/qo- 

DISCUSSION AND CONCLUSION 

The violation of the Lorentz sjTiimetry of spacetime results in "spoiling" the ordinary light 
cone and, in particular, in the emergence of birefringence in light propagation. Thus, the 
investigation of electromagnetic wave phenomena provides an understanding of the Lorentz 
violation in the photon sector of SME. The birefringence type effects represent the observ- 
able consequences of the model and the corresponding measurements impose limits on the 
Lorentz-violating parameters in the general linear spacetime relation. It was shown previ- 



ously [5[ by using cosmological and laboratory observations that the birefringence-related 
parameters of the principal part of the constitutive tensor must be smaller than 10~^^ Aq. 
In this paper we investigated the general covariant Fresnel equation for the wave covectors 
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as well as for the ray vectors in the case of a linear spacetime relation ([3]) that connects the 
electromagnetic excitation with the field strength. Strictly speaking, the physical informa- 
tion contained in the ray vector surface is the same as that in the wave covector surface. The 
reason why we still analyzed the ray along with the wave surfaces is twofold. Firstly, we de- 
veloped the four- dimensional covariant approach to the ray surfaces in complete analogy with 
the earlier results obtained for the Fresnel wave surfaces. Such a formalism might be helpful 
in the study of the various aspects of the wave propagation in media, since the construction 
of the ray surfaces is quite a common tool in crystal optics, see Q 
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3l|. Sec- 



ondly, we demonstrated explicitly that although the ray surface is dual to the wave surface 
(which is also clearly seen from the figures above), the optical metric, which is naturally 
derived in the case of birefringence, is different for the wave and for the ray surfaces and the 
two metrics derived are not dual to each other. 

The skewon part (fTO|) i . (fT5|) of the constitutive tensor was in the center of our study. We 
introduced its general parametrization ( !55|) . This enabled us to distiguish between three dif- 
ferent types of effects: the electric and magnetic Faraday effects and the (magneto-electric) 
optical activity. The influence of a skewon field on the Fresnel wave and ray surfaces is qual- 
itatively different in these three case. This is illustrated in the figures in the corresponding 
subsections of Sec. . However, the characteristic sign of the skewon is the emergence of the 
specific holes in the Fresnel surfaces that correspond to the directions in space along which 
the wave propagation is damped out completely. This effect is in c omp lete ag reement with 



our earlier conclusion on the dissipative nature of the skewon field 13|, |22|, |23 |. 



Besides the qualitative results, one may derive numerical limits on the magnitude of the 
skewon components from the experimental search for the anisotropy of the velocity of light. 
Using the data for the gamma-ray bursts and following we then obtain the estimates 
SaSfj"- ~ aj'ah'^ < 7 X and UaU'^ < 3 x lO'^'^eg. This shows that if the skewon 

indeed spoils the light cone structure (yielding a violation of the Lorentz symmetry and 
an anisotropy of the propagation of light), its influence should be extremely small. As a 
final remark we note that the experimental laboratory data 35|] imposes much stronger 
constraints on the principal part of the constitutive tensor (see the corresponding analysis 
in {sl), and we can expect that the same will be true for the skewon too since the theoretical 
formalisms are pretty much the same in both cases. 
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